Microtraps for neutral atoms using superconducting structures in the critical state 
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Recently demonstrated superconducting atom-chips provide a platform for trapping atoms and 
coupling them to solid-state quantum systems. Controlling these devices requires a full understand- 
ing of the supercurrent distribution in the trapping structures. For type-II superconductors, this 
^\ . distribution is hysteretic in the critical state due to the partial penetration of the magnetic field in 

the thin superconducting film through pinned vortices. We report here an experimental observation 
<^^> of this memory effect. Our results are in good agreement with the predictions of the Bean model of 

O^l ■ the critical state without adjustable parameters. The memory effect allows to write and store per- 

manent currents in micron-sized superconducting structures and paves the way towards new types 
^ ' of engineered trapping potentials. 

o 

Atom-chips trap ultracold neutral atoms in the magnetic field of micron-sized current-carrying wires or ferromagnetic 
r- — I- films High field gradients are obtained at short distances from on-chip microfabricated structures. The resulting 
*^ strong confinement considerably simplifies the production of Bose-Einstein condensates (BEC) @, Moreover, 
II ■ complex potentials can be designed and exploited for the realization of atomic conveyor belts Q , beam-splitters [f| , 
" pH \ and miniaturized atomic clocks jf|. 

Recent experiments 0, ID, Q have demonstrated that the normal metal wires can be replaced by superconducting 
materials. The trapping time is expected to be extremely long close to superconductors [3, [ll|, EH, EH ; due to a 
| ^ * considerable reduction of the Johnson-Nyquist noise, which is the dominant source of losses for normal metal atom- 
chips. Moreover, the cryogenic environment of superconducting chips makes it possible to couple cold ground state 
t-H ■ or Rydberg atoms with solid-state quantum devices, such as linear coplanar microwave cavities [L"j | or Josephson 
devices (l5| . These hybrid systems will open new avenues for fundamental research and have potential applications 
in quantum information processing (lfl |. 

Cold atoms on chips are also sensitive probes of the magnetic trapping potential and, hence, of the current distri- 
bution in the wires. Tiny current fluctuations caused by edge imperfections can explain atomic cloud fragmentation 
close to the chip 11711 . Sensitive imaging of the magnetic field by a BEC has been used for probing current distributions 
in metallic films [181 ]. Current distributions in superconducting films are expected to be dramatically different from 
^\ | those in normal metals. Cano et al. have recently shown that trapping frequencies decrease significantly close to a 
■ cylindrical niobium wire due to the screening currents expelling the magnetic field from the superconductor (Meissner 
effect) 0. However, this Meissner regime does not apply to superconducting thin films, where lar ge e nough magnetic 
fields partially penetrate the film through vortices. The Bean model of the critical state 1^, 2(], 21 1 provides a the- 



oretical description of the current distribution in this mixed phase. It predicts hysteresis for the film magnetization, 
■ which cannot be described by the Meissner thermodynamic phase. 

We present here an experimental evidence of this hysteretic behavior. We bring ultracold atoms close to a super- 
conducting thin strip. The atomic cloud follows trajectories which reveal different magnetic field maps depending on 
the magnetization history of the superconductor. Our observations are in good agreement with numerical calculations 
based on the Bean model [22j ]. with no adjustable parameter. These results show that permanent currents present in 
superconducting films can be understood quantitatively and controlled at the micrometer scale, opening the way to 
new trap geometries. 

The layout of our atom chip is shown in Fig. [1^ and is described in detail in Rcf. 0). The chip is based on a 
360 fjm thick silicon substrate covered by a 500 nm thick layer of silicon oxide. The structures used for trapping 
are a U-shaped wire and a Z-shaped wire. These wires are obtained by sputtering a 1 fan niobium layer (transition 
temperature T c = 9.2 K), which is then patterned using standard optical lithography. The central regions of the U- 
and Z- wires have widths of 300 fim and 40 /im respectively. The measured critical current of the Z-wire is 1.8 A; 
that of the U-wire is larger than 5 A. Following wire patterning, a 1.5 /im planarization layer is added as a base for 
a 200 nm thick gold layer. It acts as a mirror used for the mirror magneto-optical trap (MOT), which traps and 
cools the 87 Rb atoms before loading in the chip magnetic trap. The chip is vertically mounted on the cold finger of a 
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FIG. 1: (a) Central part of the chip, showing the U- and Z-wires. (b) Cross-section of the chip along the dashed line in (a). 
The initial trap position and the two approach trajectories are also represented. The position of the cloud is measured with 
respect to the origin O at the center of the Z-wire. 



custom-made cryostat, operated at a temperature of 4.2 K. 

The sequence for atom cooling and trapping is detailed in 0, [23[ . The U-wire is used in combination with an uniform 
bias field to create the quadrupolar magnetic field required for the operation of the mirror-MOT. After trapping and 
cooling in the mirror-MOT, the atoms are further cooled using optical molasses and then optically pumped to the 
\hS\/2,F = 2,mp = 2) hyperfine sublevel of the ground state of 87 Rb. We then abruptly switch on the trapping 
potential generated by a current Iz in the Z-wire and a bias magnetic field (B x , B y , B z ). The cloud is then compressed 
and evaporatively cooled. We finally adiabatically set the Z-wire current to Iz = 1-34 A and the bias magnetic field 
to (—3.0,0,9.4) G. The w 1.4 • 10 5 atoms are then at a temperature of 30(5) fiK and at a distance of about 250 ^m 
above the Z-wire. 

The final approach to the chip is controlled by linearly sweeping, during 400 ms, B y and B z towards the final values 
(B y j,B z j) while keeping B x = —3 G constant. After each sequence, an absorption image of the atomic cloud is 
taken by using a probe beam sent along an axis in the xy-plane, at an angle of 11° with respect to the x-axis. The 
absorption image shows the cloud and its reflection in the on-chip gold mirror (sec Fig. [2]). A gaussian fit of the direct 
and reflected images provides the vertical position z c and the distance y c to the chip [23j . 

We first obtain a reference trajectory by approaching the atoms to the surface while keeping a large distance from 
superconducting wires j24j. In order to image the cloud trajectory, we repeat 61 times the full sequence varying 
linearly (B y j, B z j) from (0,9.4) G to (14.1,0) G. The cloud follows the nearly circular trajectory labeled 1 in Fig.[T|D. 
The atom-Z-wire distance remains always larger than 240 [ira, while the distance y c to the chip surface decreases. 
Fig. [2^ presents an absorption image for (B y j, B z j) — (10.1, 2.7) G. We observe the cloud (on top) and its reflection. 
The cloud has a diameter of 26(3) (im and a length along the rr-axis of 280(3) /im. This finite length accounts for the 
partial overlap of the direct and reflected images. Moreover, due to the asymmetry of the Z-wire with respect to the 
yz-planc, the cloud is slightly distorted and displaced by -380(5) /jm in the cc-direction. The distance of the atoms 
to the surface is measured to be y c = 81.7(6) (im. For y c < 50 (im the atoms are rapidly lost by adsorption on the 
surface. 

We then study the influence of the U-wire by approaching the atoms along the trajectory labeled 2 in Fig. lb 
((B Vi f, B z j) vary from (0,9.4) G to (-14.1,0) G). We compare different magnetization histories for the atom-chip. 
They correspond to cooling down the superconducting wires through the transition temperature T c with different 
applied perpendicular fields B yt . Fig. [2)3 corresponds to B yt =0 G and (B y j,B z j) =(-6.6,5.0) G. The direct image 
of the cloud (on top) is less visible than the reflected one because of diffusion of the probe beam on the chip wires. For 
this magnetization history the atoms are lost when the approach-distance y c is lower than w 140 /im. This minimum 
distance is significantly larger than that of trajectory 1. 

For Figs.[2b-d, B y ^ is equal to -3.0 G and (B yJ , B zJ ) =(-5.4,5.8) G and (-8.7,3.6) G respectively. We observe that 
the cloud splits in two parts in between Figs. 2c and 2d (note that in Fig. 2d the direct and reflected images of the 
left cloud overlap). This splitting is radically different from the early cloud disappearance observed with the other 
magnetization history (Fig.[2t>). Movies of the complete cloud trajectories are available in [2^|. These approaches to 
the chip surface clearly exhibit the influence of permanent hystcrctic currents in the U-wire on the trapping potential. 
Note that the cloud configuration is found to be quite stable, leading to reproducible images over a few days, in spite 
of repeated changes of applied currents and fields. The permanent current distribution is only reset by a transition 
to the normal state. 

The current distribution in a type-II superconductor is known to be well described by the Bean critical-state 
model (l9l . 20L 26| . It considers quantities averaged over lateral dimensions larger than the intervortex distances and 
the thickness of the sample. It thus describes correctly a situation in which the distance between the atomic cloud 
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FIG. 2: Absorption images of the atomic cloud for different approaches to the chip surface. The final values of the external 
bias field are given in the text. The color scales indicate the optical density. The dotted horizontal lines separate the direct 
(top) and reflected (bottom) images of the atomic cloud, (a) Approach along trajectory 1 in Fig. [Tp. (b-d) Approach along 
trajectory 2 for B V:t ~0 G (b) and B Vt t =-3 G (c-d). (e-h) 2D-plot of the calculated trapping potential for the same conditions 
as in Figs. (a-d). The two solid lines correspond to equipotential lines at 5 and 30 /iK respectively. The red (left) and green 
(right) rectangles above the abscissa represent the extension of the Z-wire and the U-wire respectively along the z-axis. 



and the U-wire is much larger than these length scales. The main assumption of the model for a thin film is that 
the surface current density K has a modulus always smaller than a maximum value Kc- When the current tends to 
exceed locally Kc, vortices rearrange, leading to a reconfiguration so that |K| < Kc- The Bean model relies thus 
on a single parameter, Kc, which in our case can be directly inferred from the measured Z-wire critical current: 
Kc = 45 mA/fim. 

There are analytical solutions of the model in simple geometries [27l, [28j. which do not apply to our situation. In 
order to calculate the current distribution, we use thus a numerical procedure based on a variational formulation of 
the model by Prigozhin 22j. It applies to two-dimensional geometries. We assume that the central parts of the U- 
and Z-wires extend to infinity along the ^-direction. This assumption is reasonable, since the length of the central 
parts of the wires (2 mm) is much larger the atom-chip distance. Moreover, the parts of the wires parallel to the 
z-axis (see Fig. [T^,) generate a field oriented mainly in the x-direction, which does not affect to first order the trap 
position in the j/z-planc. 

With these assumptions, the surface current density K is oriented along the x-axis. The U- and Z-wires are 
discretized in elements with a 3 /xm width. Applying this discretization to a single wire, we recover accurately the 
available analytical result [28| . The complete numerical simulation takes into account the variations of the applied 
currents in the Z- and U-wires as well as those of the bias field during the whole experimental sequence. From the 
current distributions, we compute the magnetic field in the yz-planc. Figs. 2e-h present the map of the trap potential, 
proportional to the magnetic field amplitude, for the configuration of Figs. 2a-d respectively The thick solid lines 
represent the 5 and 30 /xK energy equipotentials. The predicted centers and extensions of the clouds are in good 
agreement with the observations. In the case of Fig. 2f, a shallow minimum is located on the surface. When the 
atoms are brought slightly closer to the chip, the potential barrier between the two wells lowers and the atoms rapidly 
escape towards the surface. In the case of Fig. 2h, the second minimum is always above the chip surface. As the bias 
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FIG. 3: Cloud center coordinates (y c ,z c ) as a function of the amplitude of the external bias field for (a) approach 

trajectory 1, and (b-c) approach trajectory 2 for B y .t =0 and —3 G respectively. Full (open) symbols and dashed (solid) lines 
correspond to the coordinate y c (z c ) inferred from the experimental data and from the theoretical Bean model respectively (the 
noise in the theoretical data is of numerical origin). For experimental data, it is sometimes impossible to extract one of the two 
or both coordinates from the image. This is the case when the cloud and its reflection on the chip surface are indistinguishable 
or when the signal to noise ratio becomes too small due to low atom number. 

field is varied, the two potential wells get populated. 

For a more quantitative analysis, Fig. [3] presents the comparison of the coordinates of the cloud centers (y c , z c ) as a 
function of the final bias field absolute value |-Bj/,/| together with the numerical prediction. The approach trajectory 
starts at \B y j\ ~0 G. Fig. [3^ corresponds to the reference trajectory 1. Fig. [3b-c correspond to trajectory 2 after a 
transition of the superconducting layer in a bias field B yt =0 and -3 G respectively. The experimental data (points) 
are in good agreement with the predictions (lines), with no adjustable parameter. The comparison of Figs. 3a and 3b 
shows that the presence of the superconducting layer prevents the atoms from getting close to the surface. In Fig. 3b, 
above \B y j\=6 G, we observe that the experimental distance to the surface y c is slightly larger than the theoretical 
prediction. 

In the case of Fig. 3c, theory predicts the existence of two separate potential wells above the chip surface for all 
values of B y j. For |-By,/| < 6 G, atoms arc confined to the potential minimum which is initially above the Z-wire (open 
and full circles). The two wells nearly merge around \B y j\=6 G. After merging, both wells are populated as observed 
in Fig. 2d. The coordinates of the second well are represented by open and full squares. Note that, for \B y j\ > 7 G, 
we cannot measure y c for the left cloud as direct and reflected images then merge (see Fig. 2d). As in Fig. [3b, the 
distance to the surface y c of the second well is slightly underestimated by theory. We attribute these discrepancies to 
the assumption of infinitely long U- and Z- wires along the x-direction. For large values of B y j, we measure that the 
cloud moves along x by up to 0.5 mm and gets closer to the bends of these wires. A three-dimensional calculation of 
the current distribution and trapping potential would be better suited to this situation. 

These results show that hysteretic permanent currents are an essential feature of type-II superconducting atom- 
chips. The trapping potential is quite different from that of normal structures. The experiments are in good agreement 
with numerical calculations based on the Bean model without adjustable parameter. These phenomena can be used for 
new programmable trapping geometries. A superconducting current-carrying wire guiding atoms in a one-dimensional 
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potential well can be placed between two rows of square niobium patches, whose size and spacing are in the /jm range. 
If these patches become superconducting in a bias field B Vit , permanent currents create a modulation of the trapping 
potential. These currents can be cancelled by a selective heating of the patches. Switching is fast and requires low 
energy when the current and temperature are close to the critical values (29j . It can be triggered by a weak laser 
pulse or by a resistive element integrated in the chip under the patch. Erasing the magnetization of a selected set 
of patches realizes programmable superlattices or disordered potentials on a time scale much faster than the typical 
atomic oscillation period (~ 100 /xs). This device opens interesting avenues for the study of atomic transport and 
localization. In a more complex setting, atoms trapped in a plane parallel and close to the chip by a laser standing 
wave can be laterally confined by the potential created by an array of magnetized superconducting patches. This 
programmable confinement opens the way to experiments on two-dimensional transport in complex geometries. These 
perspectives are also relevant for hybrid quantum systems 3> IHj with atomic systems interacting with solid-state 
quantum structures. 
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